Charge Ordering and Spin Dynamics in NaX^Os 
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We report high-resolution neutron inelastic scattering experiments on the spin excitations of 
NaV205. Below T c , two branches with distinct energy gaps are identified. From the dispersion 
and intensity of the spin excitation modes, we deduce the precise zig-zag charge distribution on the 
ladder rungs and the corresponding charge order: A c w 0.6. We argue that the spin gaps observed 
in the low-T phase of this compound are primarily due to the charge transfer. 



PACS numbers: 71.45.Lr, 75.10.Jm, 75.40. Gb 

The low dimensional inorganic compound NaV2 0s un- 
dergoes a phase transition at T c = 34 K |jj associated 
with both a lattice distortion Q and the opening of an en- 
ergy gap to the lowest triplet spin excitations j|] . While 
the nature of the low-T phase in NaX^Os is not fully un- 
derstood it is clear that, unlike CuGeOs, the spin-Peierls 
model does not apply simply to this compound [Q. The 
spin gap may result from charge-order (CO) rather than 
the lattice distortion ||. Indeed, NMR measurements 
indicate two inequivalent vanadium sites below T c , while 
there exists only one site above Q . There has been no di- 
rect evidence for the connection between CO and a spin 
gap, nor to distinguish various conjectured spatial dis- 
tributions of charge In this letter, we present new 
results of neutron inelastic scattering (NIS) on the spin 
excitations in the low-T phase that can now resolve these 
issues. 

In NaV205, the vanadium ions have a formal valence 
of 4.5+. Initially, this was proposed to correspond to an 
alternation of V 4+ ions, with a spin value S — 1/2, and 
V 5+ ions with S = ||. At room temperature, NaV205 
is well described by a quarter-filled two-leg ladder system, 
with only one type of vanadium site V 4 5+ . From calcu- 
lations of electronic structure d|l(J] , the strongest orbital 
overlaps are on the ladder rungs. One expects that the 
5 = 1/2 spins are carried by the V-O-V molecular bond- 
ing orbitals, with charge fully delocalized on two sites. 
As the energy of the anti-bonding orbital is much higher, 
it can be projected out, and above T c , these spins, as they 
interact in the leg direction (|| b axis), form an effective 
uniform quantum Heisenberg spin chain with interactions 
between chains that are both weaker and frustrated. 

At low temperatures, NMR shows this can no longer 
be so. On each rung, a charge transfer A c may occur. 
Taking the average charge on vanadium sites to be 1/2, 
the charges on the two vanadium sites on a rung are de- 
fined through n± — (1± A c )/2. Two forms of CO can be 



considered g], the in-line, with the same charge trans- 
fer on each rung, and the zig-zag with alternation in the 
charge along the ladders as shown in fig. lb. Recent 
X-ray diffraction measurements jll| established that the 
lattice structure below T c consists of a succession of dis- 
torted and non-distorted ladders of vanadium ions (see 
fig. la). Neglecting inter-ladder diagonal couplings J±, 
the ladders would behave magnetically as independent 
spin chains. For one ladder (chain 2 in fig. la) dis- 
tortions in the exchange paths both within the ladder 
and via neighboring ladders result in an alternation of 
the effective exchange coupling in the b direction, J^i 
and Jb2- The ladders in which the rungs are distorted 
(chains 1 and 3), however, remain magnetically uniform 
as a mirror plane passes through each rung. A minimum 
magnetic model without CO would be a succession of al- 
ternating and uniform chains. An energy gap (expected 
to be small as it results primarily from the alternation 
in JjJ would characterize the excitation branch of the 
alternating chains, and there would be no gap for the 
uniform chains. The initial NIS Q found, however, two 
excitation branches, with the same gap at the antiferro- 
magnetic point T + = E~ ea 10 mcV. To analyse these 
results, a recent spin model jl2| used an explicit relation 
between the spin excitations and the CO. This model as- 
sumed a single gap and, as it implies zero intensity of 
one excitation branch, must be extended to explain the 
NIS data. A more precise determination of excitations in 
the low-T phase of NaV205 is therefore crucial. In the 
present work, using high-resolution, the dispersion of the 
excitations is re-explored in a wider part of the reciprocal 
space. Moreover, the evaluation of the structure factor, 
i.e., the (energy- integrated) intensity of each excitation 
mode, allows us to determine the charge transfer A c . 

The single crystal (« 8x5x2 mm 3 ) was grown by a 
flux method. The NIS measurements were performed at 
T < 4.2 K, on two thermal neutron three-axes spectrom- 
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eters - IN8 and CRG/CEA-IN22 - at the Institut Laue- 
Langevin (ILL). On IN8, vertically focusing monochro- 
mator PG(002) and Cu(lll) were used in conjunction 
with a vertically focusing analyzer PG(002) and hori- 
zontal collimations 60'-40'-60'. The final wave vector 
was kept fixed at kf — 4.1 A -1 . IN22 was operated 
at kf — 2.662 A" 1 , with a PG(002) monochromator 
and a PG(002) analyzer used in horizontal monochro- 
matic focusing condition (resulting wave vector resolu- 
tion: Sq rs 0.03 r.l.u.) with no collimation. The sam- 
ple was installed in an "orange" ILL cryostat, with the 
scattering wave- vector Q lying in the reciprocal (a*,b*) 
ladder plane. 




FIG. 1. a) Simplified representation of distorted and 
non-distorted chains, chain 1 (or 3) and 2, respectively, with 
the Jj_ bond alternation between chains 1 and 2; b) Proposed 
charge order (large and small lobes represent large and small 
average charges) leading to a spin gap. The elementary trans- 
lation (a,b) of this CO agrees with the observed periodicity 
of spin excitations. 

The two branches characterizing the low-energy exci- 
tations in NaV205 have distinct energy gaps: E+ ^ E~ . 
This important result will be established when we con- 
sider the dispersions in the transverse a direction. First, 
however, we determine the dispersion in the leg direc- 
tion (b axis). Examples of constant-energy scans ob- 
tained on IN8 as a function of Qb are displayed in fig. 
2a (wave-vector components are expressed in reciprocal 
lattice units, r.l.u.). Increasing the energy, one resolves 
the single peak seen at 10 meV into propagating modes 
(dashed lines), whose peak position is shown in fig. 2b. 
They describe the dispersion of the elementary excita- 
tions in the b direction, near the AF chain wave-vector 
component (= Qb — 0.5). They are compared to 

a dispersion law characteristic of a gapped spin chain: 

E(q b ) = ^E 2 g + (El - E 2 g ) sin 2 (27rQ 6 ) where E g and 

E m are the gap and maximum energies of the disper- 
sion, respectively. For both E g = E+, E~ (the solid 
and dashed lines, respectively), we evaluate E m — 93 ±6 
meV. Compared to the prediction for a uniform Heisen- 
berg chain, E m — nJb/2, where Jb is the exchange in the 
chain, one obtains for the low-T phase, Jb ~ 60 meV. 

Second, we consider the dispersions in the a* direction 
(i.e., along the rungs). A few examples of energy scans 



performed on IN22 at constant Q are reported in fig. 3. 
In general, two peaks are observed. At a few Q a values, 
however, an extinction occurs. This extinction may con- 
cern one of the two modes, or the two modes simultane- 
ously. The left and right panels report data obtained for 
Q^ F = 0.5 and Qb = 1 (equivalent to the zone center of 
the AF chains Qf c )- As examples, we show for Q£ F that 
the smallest energy difference between the two observed 
peaks is obtained for integer Q a values (here, Q a = 3), 
the largest one for half-integer values (Q a = 2.5) and an 
extinction of the two modes occurs at Q a « 1.75 (iden- 
tical results have been obtained for Qb = 1.5). Surpris- 
ingly, at the chain zone-center Qf c we found a small but 
non-zero intensity for the two excitation branches. The 
smallest energy difference between the two peaks is now 
obtained for half-integer values (here, Q a = 1.5). Extinc- 
tions of one of the two modes is observed at Q a = 1 (on 
the high-energy mode) and at Q a — 2 (on the low-energy 
mode). For all the spectra recorded on IN22, the back- 
ground (the dotted lines in fig. 3) was carefully deter- 
mined. Several procedures have been used involving data 
recorded at low and high temperature (i.e., above T c ). 
In fig. 3a, for instance, the open dots used to define the 
background are obtained from measurements performed 
at 40 K |l3) while in fig. lb, it is determined from Qb 
scans performed at low temperature for different energies. 
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FIG. 2. IN8 data: a) Constant-energy scans as a function 
of Qt. The full lines is the result of a fit to the data, which 
takes into account both the magnetic contributions (shown 
as the dashed lines) and the background (the dotted line); b) 
Energy dispersion in the b direction. The lines correspond to 
the two branches with gaps E~£ and E~ . 

For the analysis, we assume the two observed peaks 
to belong to two distinct contributions. Their unsym- 
metrical lineshape is characteristic of gapped excitations 
undergoing a rapid energy dispersion (as established in 
fig. 2b). In such dynamical response func- 

tion (shown by the dashed lines in fig. 3) is well-suited 
to fitting, conveniently defined with only 3 parameters: 
the peak energy E + (£*-), an intensity factor A + (A_) 
and an energy damping L fl4| . As L is mainly fixed 
by the resolution conditions, it is assumed to be the 
same for the two contributions (r w 0.4 — 0.8 meV). 
Together with the background, the agreement with the 
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experiments, shown by the solid lines, is good. The val- 
ues obtained for E± as a function of Q a are shown in 
fig. 4. We establish several new features. At both Q^ F 
(solid symbols) and Q^ c (open symbols), the transverse 
dispersion consists of two distinct excitation branches 
which never cross. This justifies our previous statement, 
namely that there are two distinct energy gaps, E£ and 
E~ . In each branch, the periodicity is 27rQ a : this is 
twice that previously determined. The two dispersions 
have the same amplitude, SJ ~ 1 meV but, remark- 
ably, the upper and lower branches are out of phase, 
and there is phase inversion between branches at Qf c 
and Q^ F ■ For each excitation branch, the corresponding 
structure factors, S b AF (Q a ) and S^(Q a ) are evaluated 
by integrating the fitted dynamical response functions 
over a wide energy range (from up to E w 300r). We 
estimate systematic error for varying the upper cut off 
alters the results by at most 5%. The resulting val- 
ues (dots and squares) are reported in figs. 5a and 
b. The sums S? F {Q a ) = S£ F {Q a ) + S£f (Q tt ) and 
Sb C (Qa) = S^{Q a ) + S£F(Q ) are shown as the stars. 
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FIG. 3. IN12 data: Constant-Q scans as function of energy, 
a) at Qb = Qi, F ; b) at Qb = Qb° ■ The two arrows point to 
the positions expected for the two missing peaks. Curves are 
described in text. 

The interpretation of these results is developed in 
three steps. First, the charge order: each spin is 
associated with an electronic wave function on the 
two sites of a rung that depends on n±. The 
structure factors for the in-line and zig-zag mod- 
els are S Qb (Q a ,uj) = cos 2 (TrQ a p) S (Q a , Qb, w) + 
A 2 sin 2 {■KQaP)S (Q a + \,Q b ,uj) and S Qb {Q a ,uj) = 
cos 2 (nQap) S (Q a ,Qb,w) + A 2 sin 2 {irQap) S(Q a + \, 
Qb + respectively, with uj = E/fi and where 

p = l/a w 0.304 (I, rung length and a, lattice param- 
eter) and S (Q a ,Qb,u) is the structure factor for spins 
localized on the center of each rung . From the ratios 
of intensities for different values of momentum, one can 
extract the charge transfer A c independent of the form 
of S. In particular, we verify that the order cannot be 
in-line but predictions agree with a zig-zag order with 
A 2 » 0.35, i.e., A c w 0.6. The agreement is particularly 
good for the sums S b AF {Q a ) and S b zc (Q a )- For the abso- 



lute intensities, we need an explicit form for S which we 
take from the strongly dimerized limit (SDL), in which 
the wave function is simply a product of singlets on the 
stronger bonds. For the in-line model, for example, the 
SDL would give zero intensity at Qf c in contradiction 
with the observation (the data in fig. 5b). The in-line 
model can be ruled out. In figs. 5a and b, the predictions 
provided by the zig-zag model (solid lines) are compared 
with the experimental total structure factors S b AF (Q a ) 
and S b zc (Q a ) (solid and open stars). In fig. 5a, the 
agreement is obtained with no adjustable parameter ex- 
cept for an overall amplitude factor. Once this factor is 
determined, the results in fig. 5b depends only on A 2 . As 
can be seen, a good agreement is obtained for A 2 « 0.35. 
The low-T phase of NaV205 is very well described by the 
zig-zag model, with a rather large charge transfer. 
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FIG. 4. Dispersion of the two excitation branches in the 
a direction at Qf F and Q§ C ■ The curves are theoretical 
predictions fitted to the data (see text). 

Second, we consider the transverse dispersions. Later, 
we explain that the gap is induced by the CO. In fact, 
there are two distinct gaps because of the structural 
distortions (implying distinct couplings Jj^ 2 and 2 
as shown in fig. lb). Due to the CO, 4 different 
interchain exchange integrals must be considered, J{, 
J' 2 i J3 and J4. The two branches (associated with the 
gaps £"+ and E~) acquire a transverse dispersion, dc- 
sc ribed at Q* F and Qf c by E^ = (E+ + E~)/2 ± 

y/ (E+ - Eg) 2 / 4 + r5J 2 sin 2 (irQ a ) and E^ 

(E+ + E-)/2 ± sj (E+ - Eg) 2 /^ + SJ 2 cos 2 (ttQ q ), re- 
spectively, with SJ = J1-J2+J3-J4 In fig. 4, these 
predictions (solid and dashed lines) are compared to the 
data. Again, a very good agreement is obtained yielding 
the following evaluation E+ = 10.1 ± 0.1, Eg = 9.1 ± 0.1 
and SJ = 1.2 ±0.1 meV. The dispersion gives directly the 
alternation in the inter-ladder diagonal bonds 5 J. If wc 
assume it is dominated by the CO, we can also estimate 
the average from SJ w J1A 2 Q giving Jl ~ 2.4 meV. 

The structure factors can be calculated for each 
branch. Within the SDL approach, one obtains the con- 
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tributions S b A ^{Q a ) and Sgg{Q a ) @ shown by the dot- 
ted and dashed lines in fig. 5. In fig. 5b, the agreement 
with the data is rather good. In particular, the extinction 
phenomenon observed for each branch is well reproduced. 
In fig. 5a, while the sum S b AF (Q a ) is well described, 
we note a discrepancy between the SDL predictions and 
the individual structure factors S£F (Q a ) and S b A _ F (Q a ). 
This difficulty could be explained as follows. As shown in 
fig. lb, two successive chains are not identical. A charge 
transfer giving different average valence on the chains will 
mix intensities at Q^ F without affecting the fluctuations 
at Qfr C in qualitative agreement with the observation 
|^7| . Experimental supports for such a charge transfer 
would be useful. Our proposition for the charge order- 
ing, i.e., the zig-zag model sketched in fig. la, and our 
estimate A r « 0.6 are based on the total intensities. 
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FIG. 5. Structure factors S^±(Q a ) and S%± {Qa) for the 
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two magnetic branches in the a direction: a) at 
Qjf c . The solid and open stars represent the sums S^ F (Q a ) 
and Si C {Q a ). The data are corrected from the V 4+ atomic 
form factor f V i+ , being all normalized at Qt = 0.5. The 
dot-dashed line gives the Q a dependence of f V i+ . The other 
curves are theoretical predictions compared to the experi- 
ments (see text). 

Finally, we consider the origin of the gap. As discussed 
in the introduction, the lattice distortion alone, in iso- 
lated ladders, cannot explain the presence of two energy 
gaps and their size. To analyze the effects of the diag- 
onal couplings J±, we refer to fig. lb. Each ladder is 
seen to be a succession of two distinct clusters (shown 
by ovals in the figure) . By exact diagonalization of each 
cluster, using the effective parameters of a t-J model jl(J 
and adding a potential imposing a charge transfer, we 
evaluated Jbi and Jbi as a function of A c and the bond 



alternation d = | J&i — Jta\ / {Jbi + Jbi)- For A c w 0.6, 
the couplings underestimate the experimental value by a 
factor of about 2, but as the parameters calculated on 
the high temperature structure Jl0| ] and as such cluster 
calculations are rather crude, the agreement is satisfac- 
tory. The value obtained for the bond alternation can be 
considered as reasonable: d s=s 0.025 — 0.030. Then, using 
the experimental value Jb ~ 60 meV, one finds an energy 
gap Eg « 6 — 8 meV. This is in a fairly good agreement 
with the experimental value E g « 10 meV. This simple 
analysis supports the view that, in NaV20s, the gaps are 
primarily due to the CO fllifl . The lattice distortion plays 
a secondary role, explaining why two distinct branches 
are observed experimentally, and their separation. In our 
picture, the magnetic anisotropies are unnecessary 0. 

In conclusion, the CO in NaV20s is quantitatively de- 
termined by the present NIS measurements ]l9|. It ex- 
plains also the energy gaps observed in the low-T phase 
of this compound. 
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